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ABSTRACT 

We present high cadence (1-10 hr~^) time-series photometry of the eruptive young variable 
star V1647 Orionis during its 2003-2004 and 2008-2009 outbursts. The 2003 hght curve was 
obtained mid-outburst at the phase of steepest luminosity increase of the system, during which 
time the accretion rate of the system was presumably continuing to increase toward its maximum 
rate. The 2009 light curve was obtained after the system luminosity had plateaued, presumably 
when the rate of accretion had also plateaued. We detect a 'flicker noise' signature in the power 
spectrum of the lightcurves, which may suggest that the stellar magnetosphere continued to 
interact with the accretion disk during each outburst event. Only the 2003 power spectrum, 
however, evinces a significant signal with a period of 0.13 d. While the 0.13 d period cannot 
be attributed to the stellar rotation period, we show that it may plausibly be due to short-lived 
radial oscillations of the star, possibly caused by the surge in the accretion rate. 
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INTRODUCTION 



On 9 February 2004. lMcNeil (|2004l ) discovered a previously unknown object about 12 arcminutes south- 
west of the M78 reflection nebula. Studies of images of the area taken prior to the event confirmed that a 
star, V1647 Orionis, had brightened significantly over the course of a few days, illuminati ng the material 
surrounding it and creating what is now known as McNeil's Nebula. iBriceno et al.l (j2004l) . through their 
long-term survey of the Ori on Nebula region, con strained the onset of the outburst to early November 2003. 
The lightcurve obtained bv lAcosta-Pulido et al. I I2Q07) shows that the object reached maximum brightness 
by the beginning of March 2004 and had faded back to its init ial state by March 2006. In examining 
photographic plates from the Asiago and Harvard Observatories, Aspin et al. ( 20061) found that this star 
had under gone a similar eruption in 1966, fading back to invisibility by the end of November 1967. In 2009, 
Aspin et al. (2009) reported that yet another outburst event began in August 2008. 



The rare nature of these events in general, and the fact that we have been able to observe several of 
them from one star in particular, give us a unique opportunity to study such stellar outbursts within the 
context of star formation. Here, we present high cadence time-series photometry of V1647 Orionis from its 
2003-2004 and 2008-2009 eruptions. 
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2. LIGHT CURVE OBSERVATIONS AND REDUCTIONS 

We observed V1647 Ori during a nine-night observing run in 2003 December with the Mosaic-1 wide- 
field imager on the WIYN 0.9-m telescope at the Kitt Peak National Observatory. This instrument consists 
of eight 2048 x 4096 pixel CCDs with a plate scale of 0'.'43 pix'^ and total field of view of 59'x59'. Four of 
the nine nights were lost to poor weather. We used the SDSS z filter at 9400A to obtain a total of 65 images 
on five nights with an average cadence of ^1 hr^^. The exposure time for all images, except those taken 
on the first night, was 120 s. On the first night, we took shorter exposures that were 60 s as well as longer 
ones of 180 s. The observations were obtained at airmasses ranging from 1.8 to 2.5. This dataset samples 
the time during which the source's brightness was steeply increasing (cf. Fig.[T]). 

We used the Y4KCam on the SMARTS 1.0-m telescope at the Cerro Tololo Inter- American Observatory 
to observe V1647 Ori on UT 2009 January 23. The 4064 x 4064 CCD has a plate scale of 0'.'289 pix^^ and a 
field of view of 20'x20'. We used the Ic filter for our 32 images taken with an average cadence of ~10 hr^^. 
The exposure time for each image was 300 s and the range in airmass was 1.2 to 2.6. If we assume a rise time 
similar to the 2003 outburst, then these data were taken soon after the object peaked in brightness (Fig. [T]). 

In order to measure how the star's brightness changed with time, we performed aperture photometry 
using standard IRAF routine^. We used an aperture radius of 6 pixels and measured the sky background 
with a 5-pixel wide annulus and a 10-pixel inner radius for the 2003 data; we used an 8-pixel inner radius 
for the 2009 data. We selected these parameters based on the average seeing of the two datasets (3.6 pixels 
in 2003 and 3.0 pixels in 2009). 

We performed differential photometry because the observing conditions were non-photometric in both 
2003 and 2009. Because most of the stars in this field are likely variable, we did not attempt to choose a 
single comparison star with which to determine the differential light curve of V1647 Ori. We expect that 
on average any variations in the field stars are uncorrelated except for effects of the instrument and sky 
conditions. Thus we selected five calibration stars (Table [IJ in the vicinity of the McNeil's object that were 
of comparable brightness to V1647 Ori with which we defined an average "reference star." We checked that 
none of these are known to be variables. We further checked that the differential light curve of each of these 
stars, determined relative to the other four calibration stars, was not variable within the photometric errors. 

The full differential light curve from 2003 December is presented in Table [5] and from 2009 January in 
Table |3l The errors on individual photometric measurements are typically 0.02 mag in the 2003 data and 
0.01 mag in the 2009 data, which include both the formal photometric errors of V1647 Ori and the error of 
the mean of the combination of the five comparison stars. 



'^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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3. RESULTS 
3.1. 2003 light curve 

3.1.1. Variability 

Our data show that V1647 Ori brightened by almost 2 mag over the course of our 9-night run in 2003 
December (Fig. [21 top). This steep rise in brightness is approximately linear (in magnitudes), but with low- 
level structure superposed on top of the linear trend. To explore the structure of these low-level brightness 
variations, we de-trended the light curve as follows. 

First we subtracted a linear fit (Fig. [21 top, dashed line), which reveals a slow brightness variation with 
a peak-to-peak amplitude of ~^0.3 mag (Fig. [2l middle). These variations are reminiscent of those commonly 
observed in classical T Tauri stars, which can arise from variations in the accretion stream or from modulation 
due to star spots at the stellar rotation period, and which of ten exhibit periodic or quasi-periodic behavior 



(Type Hp and Type II, respectively, in the nomenclature of iHerbst et al.lll994l ) on timescales of ~1~10 d. 
The V1647 Ori variations appear to modulate on a timescale of 4-5 d. We cannot establish with our data 
whether this signal is strictly periodic because our light curve spans only two cycles of such a period, and 
moreover the data gaps in the light curve leave large phase gaps when the light curve is folded on such 
a period. Furthermore, we show below that this light curve modulation cannot be the rotation period of 
the star. In the following we refer to this component of the light curve variability as a "quasi-periodic" 
modulation. We defer speculation about its possible physical significance to Sec. [31 

Next we further de-trended the light curve by subtracting a best-fit sinusoid as a simple representation 
of the quasi-periodic modulation (P=4.14 d; Fig. [21 middle, dashed curve). The resulting residual light curve 
(Fig. [21 bottom) reveals very short-timescale variations that are significantly larger than the noise in our data 
(reduced chi-square is = 9.7). The amplitude of these variations is (Tims ~ 0.05 mag. 

We conducted a few simple tests to verify that none of the periodic photometric signals discussed here 
and in what follows correlate with seeing variations. We performed aperture photometry on a portion of 
the nebula itself using, as before, a 6 pixel aperture radius; we separately performed aperture photometry 
on V1647 Ori using a larger (9 pixel) aperture radius. We executed a periodogram analysis (as below) on 
the resultant light curves and also on the seeing variations, which were obtained by measuring the average 
FWHM of our calibration stars in each image. We found no significant differences in the target light curve, 
and no evidence for significant periodicities in either the nebular light curve or in the seeing variations (and 
in particular not at the periods reported below). We conclude that changes in seeing are not driving the 
periodic photometric variations we report in this work. 



3.1.2. Periodogram analysis 

To examine the light-curve variations of V1647 Ori in detail, we subjected the 2003 light curve (Fig. [2}) 
to a standard Lomb-Scargle power spectrum analysis. The Lomb-Scargle periodogram is well-suited to un- 
evenly sampled data such as ours. It moreover possesses well characterized statistical properties that permit 



quan titative assessment of the statistical significance of any periodic behavior in the data (see [Press et al 



I992L and references therein). We will exploit these statistical properties below. 

Fig. [31 shows the power spectrum of the non-detrended light curve (Fig. [21 top) over the frequency range 
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0.1-18 d""'^. The low frequency cutoff corresponds to 1/T where T is the total timespan of the data while 
the high frequency cutoff corresponds to half the sampling frequency (i.e., the Nyquist limit). Overall the 
power spectrum rises toward smaller frequencies, with a slope that closely approximates a l/-\/7 dependence 
(represented by dashed/dotted lines in the figure). In addition, the power spectrum exhibits several peaks on 
top of the l/v7 slope. The broad peak at / ^ 0.2 d^^ corresponds to the slow, quasi-periodic modulation 
discussed in Sec. 13.1.1] (see also Fig. [2l middle). The two peaks near / = 0.8 d~^ and / = 1.2 d""'^ are aliases 
of the / ^ 0.2 d~^ modulation beating against the diurnal data gaps (/ = 1 d^^). The peak near / = 8 d~^ 
is due to the short-timescale variability in the detrended light curve (Fig. [2] bottom). 

The power spectrum of the detrended light curve is shown in Fig. |4] (top, black curve). Not surprisingly, 
nearly all of the power at low frequencies has been eliminated by the de-trending of the linear rise and 
of the slowly varying quasi-periodic modulation. The power spectrum shows a strong peak at 7.7 d~^, 
corresponding to a period of 0.13 d, and several other strong features at nearby frequencies, which have 
power levels corresponding to a statistical confidence of 90% or higher (see below). When we filter out the 
0.13-d period peak by subtracting the best fitting sinusoid from the light curve, all of the other statistically 
significant peaks in the periodogram are also removed (Fig. 21 top, blue curve), showing them to be aliases 
and beats of the 0.13 d period. The light curve (Fig. [21 bottom) is shown folded on this period in Fig. [H 
[bottom), with the best-fitting sinusoid overlaid in blue. The amplitude of this sinusoid is 0.051 mag. 

To empirically determine the false-alarm probability (FAF) corresponding to different power lev els in 



the periodogram, w e used a Monte Carlo boo tstrapping technique as described in iPress et al.l (|l992l ) and 



implemented in e.g. Istassun et aD |l999 . I2OO4). We generated 10,000 artificial light curves by shuffling the 



actual measurements in temporal order and sampling at the same timestamps as the actual data. In this 
way, the artificial light curves retain both the noise properties and the time windowing of the real data. For 
each of the 10,000 artificial light curves, we calculated a power spectrum as for the real data and recorded 
the power level of the strongest peak in each. The resulting distribution of these 10,000 maximum peak 
heights gives directly the probability of a given peak height occurring by chance. 

The distribution of maximum peak heights is shown in Fig.jSlfor 10,000 artificial light curves simulating 
the detrended 2003 light curve (Fig. [2l bottom). From the figure, peak heights with power levels above ~7 
occur in fewer than 10% of the simulated power spectra, whereas peaks with power levels above ~11 occur in 
fewer than 0.1% of the simulated power spectra; these then define the 10% and 0.1% FAF levels, respectively 
(or equivalently, the 90% and 99.9% confidence levels). These confidence levels are represented by horizontal 
dotted lines in the observed power spectrum (Fig. [H top) . 



As described by iPress et al.l (|1992I ) , the FAP for a given peak height in a Lomb-Scargle periodogram is 
expected to follow an analytic relationship (cf. their Eq. 13.8.7) shown in Fig. [Has a dotted curve. Evidently, 
our data and its associated power spectrum closely follow the expected statistical behavior. Fig.jSlshows the 
FAP calculation for the 0.13 d period. The observed peak height (shown as vertical dashed line) has a FAP 
of 1.3 X 10^^, and is therefore very highly statistically significant. 

We checked that the 0.13 d period and its statistical significance are not dependent on the details of the 
sinusoidal de-trending that we performed in Sec. 13.1.11 As a simple alternative to the sinusoidal detrending 
(see Fig. 121 middle), we instead simply shifted all of the data points from a given night such that the mean 
differential z magnitude for the observations made on each night was 0.0. We then performed the same 
periodogram analysis as above. We recovered the same 0.13 d period as before with a FAP of 3.6 x 10"^, 
again highly statistically significant. 
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3.1.3. Stellar Rotation 



Since periodic signals in young, low-mass stars are often associated with rotation, we explored the 
possibility that we have detected stellar rotation in our light curve . The rotation per iods of low-mass pre- 
main-sequence stars are most typically in the range ~2-10 d fe.g.. iHerbst et al.ll2002l). though so me young 



low-mass stars have been observed to rotate with periods as short as ~ 0.1 d (jStassun et al 



Adopting values for the visual extinction (19±2 mag), K-band veiling (1.5±0.2 mag), spectra l type (MOzfc 



2 sub classes), and K-band apparent magnitude (9.9 mag in February 2007) as determined bv lAspin et al 
(120081) ■ a bolometric correction of BCv = —1.3 and a, {V — K) color of 3.7 as appropriate for an MO spectral 
type (Kenvon fc Hart mannlll995| ). the equations for ab solute K magnitude and bolometric luminosity from 



Greene fc Lada (,1997.) . and a distance of 426 ± 20 pc (jMenten et al 
20071 ). we calculate a stellar radius of 4.2+};° R©- 



2007 



Kraus et al. 



2007 



Hirota et al 



Aspin et al.l (|2009l ) observed line broadening in the spectrum of V1647 Ori of 120 km s^^. While the 
spectra of EXor and FUo r eruptive variables can be dominated by the hot "atmosphere" of the inner accretion 
disk during outburst, the I Aspin et al. I (|2009l) observations were obtained during quiescence of the system; we 
therefore presume that the observed spectral broadening is stellar in origin. This then provides a lower limit 
on the stellar rotational velocity of wsini=120 km s"^ . We can also assume an upper limit from break-up 
considerations of Urot~190 km s^^. 

With these upper and lower bounds on the rotational velocity of V1647 Ori, if we simultaneously push 
all other measured parameters to their 1 cr limits so as to produce the shortest and longest possible rotation 
periods, we find with greater than 99.9% confidence that the stellar rotation period lies between 0.6 and 2.7 
days. To test the robustness of our calculations, we also calcula ted the stellar r a dius u sing the quiescent 
V-band extinction and K-band apparent magnitude obtained by [Abraham et al.l (|2004l ) prior to the 2004 
outburst (~ 13 and 10.3 magnitudes, res pectively): from this, and assuming an inclination angle of 61° ± 14° 
as found bv lAcosta-Pulido et al.l (|2007l ). we obtain a most likely rotation period of ~1 d, consistent with 
the above range. If this is indeed the rotational period of the star, the diurnal gaps in our lightcurve 
would preclude its detection. In any case, we conclude that the 0.13 d period, and the ~4 d quasi-periodic 
modulation, do not correspond to the rotation period of V1647 Ori. 



3.2. 2009 light curve 

Figure [S] displays our 2009 light curve. The dataset was obtained at very high cadence 10 hr^^) over 
a timespan of several hours on a single night and, thus, is sensitive only to high frequency variations. A 
simple test shows the light curve to be variable, with reduced = 10.0. However, applying the same 
periodicity analysis as above, we find no significant periodicity in this light curve across the frequency range 
to which we are sensitive, from 7 d^^ to 135 d^^; the strongest peak we obtain has a FAF of 41%. 

The full duration of the 2009 light curve is (coincidentally) 0.13 d (7.7 d^^), and thus can be used to 
check for the presence of a 0.13 d period during the 2009 observations. As shown in Fig. |6l a 0.13 d period 
such as that observed in the 2003 light curve (Fig.|4}D) is not present in the 2009 light curve. As an additional 
check, we injected a 0.13 d sinusoid of varying amplitude into these data and found that we were able to 
induce a significant peak in the power spectrum (FAP < 1%) only if the amplitude of the sinusoid is larger 
than 0.027 mag. We thus conclude that there is no evidence for an 0.13 d period in the 2009 data with 
an amplitude greater than 0.027 mag, and we can definitively rule out an 0.13 d signal with an amplitude 



- 6 - 



of 0.05 mag as seen in the 2003-2004 outburst light curve (Fig. [2]). To summarize, we find that this 2009 
hghtcurve is variable, but no periodic phenomenon with a frequency in the range to which we are sensitive 
(7 d^^ to 135 d^^) is driving this variability. 

In Fig. [3] we show the power spectrum of the 2009 light curve data together with that from the 2003 data. 
The two datasets sample a mostly disjoint range of temporal frequencies. However, the two power spectra 
together are consistent with a single 1/%// behavior for the power spectrum as a whole. We suggest one 
possible interpretation for this l/\ff behavior below (Sec. 14. 2|) . The only clearly evident difference between 
the 2003 and 2009 power spectra is in the overlap region near / = 10 d~^; as discussed above the 2003 light 
curve exhibits a strong 0.13 d period whereas the 2009 light curve does not. 



4. DISCUSSION 

We have found strong evidence to suggest that V1647 Ori exhibited a highly significant 0.13 d periodicity 
in brightness during the rapid brightening phase of its 2003 outburst. The amplitude of this periodic variation 
was ^0.05 mag. The presence of this feature in the 2003 data, obtained while the object was in the brightening 
phase, and its absence in the 2009 data, obtained after V1647 Ori had peaked in brightness, suggest that 
this phenomenon, whatever its cause, is associated with the unstable period of time when the brightness of 
V1647 Ori was most rapidly increasing. We have also found evidence for a l/-\/7 slope in the power spectrum 
of V1647 Ori over a large range of temporal frequencies, 0.1 < / < 135 d^^. 

In this section, we consider whether the 0.13 d period may be ascribed to short-lived stellar pulsations, 
perhaps triggered by the high accretion event. Next, we discuss "flickering" in the light curve, evidenced by 
the 1 / ^f/J slope in the power spectrum, in the context of a magnetically channeled accretion flow. Finally, we 
consider whether oscillations originating in the accretion disk, similar to those observed during cataclysmic 
variable (CV) star outbursts, could be observed in young star outbursts, including FUor and EXor events, 
and we speculate that the ^4 d quasi-periodic modulation observed in the 2003 light curve could correspond 
to such an oscillation. The consideration of these mechanisms here is speculative; our aim is to examine 
whether these explanations may be plausible, but we cannot yet establish that these are definitive driving 
mechanisms for the observed variability. 



4.1. Stellar pulsations 



Since it is inconsistent with the likely rotation period of V1647 Ori (Sec. I3.1.3p . we explore the pos- 
sibility that the 0.13 d period mig ht be a manifesta tion of pulsation. Given its effective temperature and 
estima ted mass of 0.8 ± 0.2 M(7) (Aspin et al. 20081) . and based on comparison with stellar evolutionary 
tracks (jP'Antona fc Mazzitellilll994[ ). we find that V1647 Ori may lie just within the theoretically predicted 
deut erium instatj ility strip. The expected fundamental mode pulsation period would be approximately 0.5 d 
re.g.. lTomalll97i . somewhat longer than the period we detect. In addition, the 0.13 d period does not appear 
in the 2009 data. It is extremely unlikely that VI 647 Ori transitioned from being within the instability strip 
to being outside the instability strip between our 2003 and 2009 observations. 



We next consider whether the dramatic increase in the accretion rate could have induced short term 
radial oscillations of the stellar surface. From the idealized homogeneous compressible model presented by 
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Cro3 (|l980l) . the period of oscillation, for purely radial pulsation, varies according to 



P'^GM 



-4 + ri(2n2 + 5n + 3) 



(1) 



where i? is the star's radius, P is its pulsational period, G is the gravitational constant, M is the stellar 
mass, Fi the adiabatic exponent, and n the pulsational mo de. For our purposes here, this idealized model is 
not significantly different from more sophisticated models (jTassoul and Tassoullll968l ). Adopting the stellar 
parameters of V1647 Ori (see Sec. 13.1. 3p and assuming as one extreme assumption that no part of the star 
is ionized (Fi = |), we find that the 0.13 d period could correspond to pulsational modes ranging from 2.8 
to 4.7. In order to determine the effect of ionization zones on our results, we looked at the extreme case of 
a fully ionized gas (Fi = |); in this case we find pulsational modes between 3.3 and 5.4. Hence, any errors 
in not taking ionization zones into account do not significantly change our results. We therefore infer that 
the 0.13 d period most likely corresponds to a radial pulsation in the 4*'' oscillation mode, but the 3'"'^ and 
5*'' modes are also permitted within the observational uncertainties in the stellar properties of V1647 Ori. 

Pulsation solely in a h igher qyerton e radial mode, although rare, is strongly dependent on the location 
of the driving mechanism ([Kurtz 2006). Examples of si milar short-term oscillation-producing phenomena 
include recent solar observations ( Karoff fc Kieldsenll2008h that revealed that energetic magnetic reconnection 
and coronal mass ejection events can induce short-lived high-freque ncy oscillations of the s olar surface. In the 
Sun, these are confined to the vicinity of the triggering event (e.g., Kienreich et al. 20091) . The triggering of 
radial oscillations by the sudden onset of highly energetic accretion on V1647 Ori thus provides one possible 
explanation for the observed variability. 



4.2. Flickering 



Flickering is def ined as random, small amplitude brightness variations recurring on dynamical timescales 
(jKenvon et al.l 120001 ). Sometimes interpreted as an observable consequence of an inhomogeneous accretion 
flow, flickering has been observed in CV stars and, as it is associated with ac cretio n, could be obs e rvable 
during similar outburst events in young stars. Indeed, Kenvon et al. ( 2000l) and Rucinski et al. ( 20081 ) 
observed flickering in the lightcurves of FU Ori and TW Hya, respectively, the latter through observation 
of a l/\/7 slope in the power spectrum. The combination of a high observing cadence and a relatively long 
time baseline means that our data are sensitive to a large range of frequencies, and, as such, any signs of 
this phenomenon should be readily apparent in an analysis of the power spectrum of our data. 

As discussed in Sec. [31 the overall power spectrum of the observed phases of the 2003-2004 and 2008- 
2009 outbursts of VI 647 Ori follow a \/\/J trend (Fig. [3]), which suggests that flickering is one possible 
origin of the random variability components of the light curve. 

The nature of such "flicker noise" is not entirely understood; however, it appears to be linked to situations 
in which a flow (in this case the accretion of material f rom disk to star) is tunneled or is in some way forced to 
pass through a physically confined region |Presslll978l ). Indeed, accretion in young low-mass stars is typically 
envisioned to occur via magnetospheric tunneling of disk material along stellar field lines that thread the 
inner accretion disk (e.g. lShu et al.lll994 ). Also, the flickering observed in th e CV system T CrB, fo r example, 
has been attributed to turbulence in the inner regions of the accretion disk ([Pobrotka et al.ll2010l ). Thus the 
observation of flickering in V1647 Ori might suggest an interaction between the stellar magnetosphere and 
the inner regions of the accretion dis k during outbu rst, with material continuing to accrete onto the star 
along stellar magnetic field lines fe.g. IShu et al.lll994[ ). 
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4.3. Dwarf Nova-Like Oscillations in the Keplerian Inner Accretion Disk? 



We have observed a modulation with a timescale of ^4 d in the 2003 light curve of V1647 Ori (Fig. [2]). 
As discussed, our data do not permit us to determine whether this modulation is strictly periodic, or indeed 
whether it persists for more than ^2 cycles. It is nonetheless a potentially interesting quasi-period that 
could arise in a number of different ways. 

Dwarf-nova oscillations (DNOs) are quasi-periodic brightness variations typically observed in cata- 
clysmic variable st ar outbursts. These phenomena are associated with accretion, and they may persist 
during quiescence (jPretorius et al.ll2Q06l ). The oscillations are not observed in all CV outbursts, and the 
oscillation frequencies may change during a high- accretion event. Because of their relatively short periods, 
DNOs are believed to originate at the inner edges of Keplerian accreti on disks. They have only rarely been 
observed during the rise of a CV outburst ( Warner fc Pretoriui 20081) . but their behavior is thought to be 
well described by a low-inertia magnetic acc retor model in which ac cretion induces variations in the angu- 
lar velocity of an equatorial accretion belt (| Warner fc Woudti 120021 ). In addition, so-called quasi-periodic 
oscillations (QPOs), which are longer term, less coherent oscillations, have also been observed during some 
CV outbursts. There is no generally accepted model fo r what causes the QPO brightness modulations, but 
they are thought to be accretion disk phenomena (e.g.. IWoudt fc Warner 2002 ) . In any e vent , an empirical 
relationship links DNOs to QPOs such that Pqpq « 15 x Pdno | Warner fc Woudti liooj) . 



The low-inertia magnetic accretor model predicts that the DNO oscillation quasi-period (Pdno) increases 
as the accretion rate (M) decreases. It also predicts that Pdno corresponds to the Keplerian orbital period 
of the inner edge of the ac cretion disk. If we assume a Pdno~-^ relation such as th at found empirically fo r 
the dwarf nova SS Cyg bv iMauchel (|l997l ) and use the M found for V1647 Ori bv iMuzerolle et al.l (|2005l ). 
we find Pdno ~ 5.7 d. We note that the value of M might differ for our data, bringing Pdno closer to our 
observed quasi-period of 4.14 d. For example, using the Pdno~-^ relation, we calculate that M ~ lO""* 
M0 yr~^ would yield a period of 4.2 d. In any event, if this model applies to the V1647 Ori system, then, 
taking the mass of the central star to be ~ 0.8 M©, we calculate that the inner edge of the accretion disk 
surrounding V1647 Ori was located ~2.5 stellar radii from the star prior to the peak of the outburst. This 
model could thus suggest that the stellar magnetosphere was still capable of holding off the inner disk, at 
least during the early stages of the 2003-2004 outburst. 

Additionally, using the relation of Pqpo ~ 15 x Pdno, we would expect to observe a QPO with a period 
of ~60 d, which is very close to the 56 d period found by lAcosta-Pulido et al.l (|2007l) who attributed it to 
dense circumstellar clumps orbiting the star. However, given the non-uniform sampling of our data and the 
fact that these oscillations are observed over only ~2 cycles, the reality of this period, and the link between 
these periodicities and those observed in white dwarf stars, is unclear. Nevertheless, they are potentially 
intriguing. Observations probing similar timescales during FUor/EXor outbursts would allow us to further 
investigate the presence of such phenomena in young star outbursts. 



5. SUMMARY 

In this work, we present high cadence time-series photometry of VI 647 Orionis during its two most recent 
outburst events. The overall power spectrum of our 2003 and 2009 datasets displays a l/-\/7 'flicker noise' 
spectrum. The detection of 'flicker noise' in the power spectra of our light curves suggests that accretion 
continues to be mediated by the stellar magnetosphere during the observed phases of the 2003-2004 and 
2008-2009 outbursts. This picture is bolstered by the observation of a quasi-periodic modulation in the 2003 
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light curve with a timescale of ^4 d, perhaps arising from a CV-hkc quasi-periodic oscillation of the inner 
accretion disk edge at a height of 2 3 stcUar radii from the stellar surface. 

Our Fourier analysis of our 2003 detrended lightcurve, obtained mid-outburst, yields a periodic variation 
on a timescale of 0.13 d that persists in spite of the dramatic rise in brightness caused by the outburst. This 
period, detected at very high statistical significance, is not attributable to the expected rotation period 
calculated from other measured properties of the star. This 0.13 d period is absent from our 2009 light 
curve. 

The 0.13 d period is very coherent in the 2003 dataset, and it is therefore likely to have been a truly 
periodic phenomenon during those observations. Given that we do not detect this period in our 2009 
light curve, obtained post-outburst, we conclude that it is probably an accretion-induced process associated 
with the epoch when the brightness of the object is increasing (and so perhaps when the accretion rate is 
increasing), one likely candidate being short-lived radial pulsations of the star. 
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Fig. 1. — Our datasets in context. The 2003 data were taken in 2003 December, during the phase of steepest 
increase in brightness. No outburst hghtcurve has yet been pubhshed for the 2008-2009 event, but if we 
assume a rise-time similar to the 2003-2004 outburst, our 200 9 data would sample the p hase when the 
brightness of V1647 Ori had plateaued, as shown. Adapted from lAcosta-Pulido et al.l (|2007l) . 
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Table 1. J2000 coordinates of stars used for differential photometry 



RA" 


Dec*" 


05:46:22 


-00:03:37 


05:46:29 


-00:03:47 


05:46:31 


-00:04:21 


05:46:28 


-00:09:59 


05:46:26 


-00:10:25 



^in hh:mm:ss 
'^in dd:inm:ss 
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1792 1794 1796 

JD (+2451187) 



1798 



1800 



Fig. 2. — (Top): Differential light curve of V1647 Ori obtained mid-outburst in 2003 December. The system 
brightened by '-^ 2 mag over the course of 9 nights. For reference, the dotted line represents a linear trend fit 
to the data. The scale o f the time axis has b een chosen to allow direct comparison with the long-term light 
curve data presented in iBriceno et al.l (|2004 cf. their Fig. 3). (Middle): Same light curve de-trended with 
the linear trend from top, revealing a quasi-periodic modulation with a timescale of ^4 d. For reference, the 
dotted curve represents a sinusoidal trend fit to the data (P=4.14 d). (Bottom): Same light curve further 
de-trended with the sinusoidal trend from middle, revealing short-timescale variations with amplitude ~0.05 
mag. 
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Fig. 3. — Power spectrum resulting from a Lomb-Scargle periodogram analysis of the non-detrended light 
curve of V1647 Ori, on a log-log scale. The frequency range 0.1 < / < 18 d~^ (black) is provided by the 
2003 light curve data (Fig. [2l top), whereas the frequency range 7 < / < 135 d^^ (red) is provided by the 
high-cadence 2009 light curve data (Fig. [5]). Following the usual definition of the Lomb-Scargle periodogram, 
the ordinate gives the spectral power normalized by the variance of the data. The power spectrum overall 
follows a slope of (dashed/dotted lines), and in addition exhibits peaks near frequencies of ~ 0.2, 0.8, 

1.2, and 8 d^^ (see the text for discussion of the meaning of these peaks). 
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Fig. 4. — Top: Power spectrum of the detrended light curve from Fig. [Jt (black). The highest peak is at 
7.7 d~^ (a period of 0.13 d). Dotted lines represent the peak heights corresponding to confidence levels of 
90% and 99.9%; a peak above the 99.9% confidence line, for example, would have a false-alarm probability 
(FAF) lower than 0.1%. The signal at 0.13 d, is very highly statistically significant, with a FAP of 1.3 x 10^^. 
(See the text and Fig. [5]for FAP details.) The blue curve shows the result of removing the 0.13 d period; no 
significant periods remain. Bottom: Light curve data are phased on the 0.13 d period with best-fit sinusoid 
in blue. The amplitude of the sinusoid is 0.051 mag. 
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Fig. 5. — The false alarm probability (FAP) of the 0.13 d period as a function of peak height as determined 
from Monte Carlo simulations. The solid line represents the probability determined from our d ata, and the 



dotted curve, which nearly overlays the solid line, shows the analytical probability described bv lPress et al 



(jl992l ). The vertical dashed line shows the peak height we actually observe in Fig.|4^. The FAP of the 0.13 d 
period is 1.3 x 10^^. 



- 17 - 




854.64 854.66 854.68 854.70 854.72 854.74 854.76 854.78 
JD (2454000+) 



Fig. 6. — Lightcurve of V1647 Ori during 2008-2009 outburst, spanning 0.13 d. These data were taken after 
the object had reached maximum brightness. No obvious periodic brightness variations are evident. The 
best-fit 0.13 d sinusoid observed in the 2003 outburst data (Fig. |4}3) is overplotted here for comparison. 



- 18 - 



Table 2. Differential photometric measurements of V1647 Ori in 2003 December 



2978 


.7040 





.000 


0, 


.046 


2978 


.7100 





.020 


0, 


.032 


2978 


.7410 


-0 


.034 


0, 


.051 


2978 


.7440 


-0 


.026 


0, 


.032 


2978 


.7560 


-0 


.020 


0, 


.048 


2978 


.7590 


-0 


.073 





.028 


2978 


.7710 


-0 


.005 





.047 


2978, 


.7740 





.014 


0, 


.029 


2978, 


.8000 


0, 


.095 


0, 


.046 


2978, 


.8040 


0, 


.056 


0, 


.029 


2978, 


.8300 


0, 


.194 


0, 


.049 


2978, 


.8330 


0, 


.082 


0, 


.028 


2978, 


.8580 


-0 


.012 





.041 


2978, 


.8620 


0, 


.045 


0, 


.028 


2978, 


.8890 


-0 


.021 





.044 


2978, 


.8930 


-0 


.026 


0, 


.025 


2978, 


.9200 


0, 


.150 


0, 


.050 


2978, 


.9230 


0, 


.101 


0, 


.028 


2978, 


.9510 


0, 


.182 


0, 


.053 


2978, 


.9540 


0, 


.021 


0, 


.027 


2978, 


.9790 


0, 


.021 


0, 


.053 


2978, 


.9810 


-0 


.019 





.030 


2979 


.0060 


-0, 


.089 





.053 


2979 


.0090 


-0, 


.028 





.034 


2982 


.7000 


-0 


.661 





.030 


2982 


.7320 


-0 


.725 





.024 


2982, 


.7610 


-0 


.648 


0, 


.028 


2982 


.7900 


-0 


.782 





.022 


2982 


.8190 


-0 


.863 





.022 


2982 


.8600 


-0 


.769 





.020 


2982 


.8890 


-0 


.673 





.023 


2982 


.9850 


-0 


.687 


0, 


.026 


2983 


.6870 


-0 


.690 


0, 


.022 


2983 


.7150 


-0, 


.729 


0, 


.021 


2983 


.7440 


-0 


.723 





.019 


2983 


.7720 


-0 


.719 





.019 


2983 


.8000 


-0 


.705 





.018 


2983, 


.8290 


-0, 


.713 


0, 


.018 


2983, 


.9340 


-0 


.656 


0, 


.021 


2983, 


.9650 


-0 


.790 





.020 


2983 


.9940 


-0 


.833 


0, 


.021 


2984 


.6570 


-0 


.998 


0, 


.017 


2984 


.6710 


-0 


.963 





.018 


2984 


.7200 


-0 


.939 





.016 


2984 


.7490 


-1 


.011 





.015 


2984 


.7780 


-1 


.063 





.015 


2984, 


.8060 


-1 


.045 


0, 


.014 


2984, 


.8350 


-0, 


.981 


0, 


.015 


2984, 


.8730 


-1 


.081 





.014 


2984, 


.9020 


-1 


.149 


0, 


.015 


2984, 


.9300 


-1 


.056 





.016 


2984, 


.9600 


-1 


.068 





.016 


2984, 


.9910 


-1 


.033 





.017 


2986 


.6640 


-1 


.569 





.011 


2986 


.6920 


-1 


.636 





.009 
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Table 2 — Continued 



JD* Az cr, 



2986 


.7210 


-1 


.668 


0. 


,010 


2986 


.7490 


-1 


.695 


0. 


,009 


2986 


.7780 


-1 


.619 


0. 


,009 


2986 


.8060 


-1 


.589 


0. 


,010 


2986 


.8340 


-1 


.674 


0. 


,008 


2986 


.8620 


-1 


.661 


0. 


,009 


2986 


.8910 


-1 


.592 


0. 


Oil 


2986 


.9240 


-1 


.570 


0. 


,012 


2986 


.9.520 


-1 


.599 


0. 


,012 


2986, 


.9800 


-1, 


.635 


0. 


,013 



"Julian Date (+2450000). 



Table 3. Differential photometric measurements of V1647 Ori in 2009 January 



JD* AI a, 



4854. 


.6518 


0. 


.023 


0, 


.008 


4854. 


.6576 


0. 


.000 


0, 


.006 


4854. 


.6628 


0. 


.017 


0, 


.006 


4854. 


.6665 


-0. 


.019 


0, 


.007 


4854. 


.6702 


0. 


.015 


0, 


.008 


4854. 


.6739 


-0. 


.010 


0, 


.008 


4854. 


.6776 


-0. 


,000 


0, 


.007 


4854. 


,6854 


0. 


.000 


0, 


.009 


4854. 


.6891 


0. 


.037 


0, 


.009 


4854. 


.6927 


-0. 


.009 





.008 


4854. 


,6964 


-0. 


,003 





.008 


4854. 


,7003 


-0. 


,030 





.010 


4854. 


,7042 


0. 


.019 





.007 


4854. 


,7079 


-0. 


.006 


0, 


.009 


4854. 


,7116 


0, 


,061 


0, 


.009 


4854. 


,7153 


0, 


,019 


0, 


.009 


4854. 


,7190 


0, 


,045 


0, 


.010 


4854. 


,7229 


-0. 


Oil 





.009 


4854. 


,7266 


-0. 


,047 





.008 


4854. 


7302 


-0. 


,022 





.010 


4854. 


7339 


0, 


,046 


0, 


.008 


4854. 


7376 


-0. 


,063 





.009 


4854. 


7414 


0. 


,022 


0, 


.012 


4854. 


,7451 


-0. 


,031 


0, 


.008 


4854. 


,7487 


0. 


,013 


0, 


.012 


4854. 


,7524 


0. 


,043 


0, 


.011 


4854. 


7561 


-0. 


,030 





.011 


4854. 


7601 


0, 


,002 


0, 


.010 


4854. 


7637 


0, 


,008 


0, 


.011 


4854. 


7674 


0, 


,002 


0, 


.012 


4854. 


,7711 


-0. 


,007 


0, 


.011 


4854. 


,7748 


-0. 


,040 


0, 


.012 



"Julian Date (+2450000). 



